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Singularities in (3 + l)-dimensional Horava-Lifshitz (HL) theory of gravity are studied. These 
singularities can be divided into scalar, non-scalar curvature, and coordinate singularities. Because of 
the foliation-preserving diffeomorphisms of the theory, the number of scalars that can be constructed 
from the extrinsic curvature tensor Kij, the 3-dimensional Riemann tensor and their derivatives is 
much large than that constructed from the 4-dimesnional Riemann tensor and its derivatives in 
general relativity (GR). As a result, even for the same spacetime, it may be singular in the HL 
theory but not in GR. Two representative families of solutions with projectability condition are 
studied, one is the (anti-) de Sitter Schwarzschild solutions, and the other is the Lu-Mei-Pope 
(LMP) solutions written in a form satisfying the projectability condition - the generalized LMP 
solutions. The (anti-) de Sitter Schwarzschild solutions are vacuum solutions of both HL theory and 
GR, while the LMP solutions with projectability condition satisfy the HL equations coupled with an 
anisotropic fluid with heat flow. It is found that the scalars K and KijK 13 are singular only at the 
center for the de Sitter Schwarzschild solution, but singular at both the center and r = (3Af/| Al) 1 ^ 3 
for the anti-de Sitter Schwarzschild solution. The singularity at r = (3M/IAQ 1 ' 3 is absent in GR. 
In addition, all the generalized LMP solutions have two scalar curvature singularities, located at 
either r — and r = r s > 0, or r = n and r = r-z with r2 > ri > 0, or r = r 3 > and r = oo, 
depending on the choice of the free parameter A. 

PACS numbers: 04.60.-m; 98.80.Cq; 98.80.-k; 98.80.Bp 



I. INTRODUCTION 

There has been considerable interest recently on a the- 
ory of quantum gravity proposed by Horava moti- 
vated by Lifshitz theory in solid state physics [2], for 
which the theory is usually referred to as the Horava- 
Lifshitz (HL) theory. The HL theory is based on the 
perspective that Lorentz symmetry should appear as an 
emergent symmetry at long distances, but can be fun- 
damentally absent at high energies Q. With such a 
perspective, Horava considered systems whose scaling 
at short distances exhibits a strong anisotropy between 
space and time, 

x->fx, t^£ z t. (1.1) 

In (3 + l)-dimensional spacetimes, in order for the the- 
ory to be power-counting renormalizablc, it needs z > 
3. At low energies, the theory is expected to flow to 
z = 1 , whereby the Lorentz invariance is "accidentally re- 
stored." So, the HL theory is non-relativistic, ultra-violet 
(UV) complete, explicitly breaks Lorentz invariance at 
short distances, but is expected to reduce to general rel- 
ativity (GR) in the infrared (IR) limit. 

The effective speed of light in this theory diverges in 
the UV regime, which could potentially resolve the hori- 
zon problem without invoking inflation The spatial 



curvature is enhanced by higher-order curvature terms 
P-0| , and this opens a new approach to the flatness prob- 
lem and to a bouncing universe [ESQ- In addition, in 
the super-horizon region scale-invariant curvature per- 
turbations can be produced without inflation [Tol . fTT| . 
and the perturbations become adiabatic during slow-roll 
inflation driven by a single scalar field and the comov- 
ing curvature perturbation is constant (Tlj . Due to all 
these remarkable features, the theory has attracted lot of 
attention lately pHbj. 

To formulate the theory, Horava assumed two condi- 
tions - detailed balance and projectability (He also con- 
sidered the case where the detailed balance condition 
was softly broken) [l[. The detailed balance condition 
restricts the form of a general potential in a (D + 1)- 
dimensional Lorentz action to a specific form that can 
be expressed in terms of a D-dimensional action of a rel- 
ativistic theory with Euclidean signature, whereby the 
number of independent-couplings is considerably limited. 
The projectability condition, on the other hand, orig- 
inates from the fundamental symmetry of the theory - 
the foliation-preserving diffeomorphisms of the Arnowitt- 
Deser-Misner (ADM) form, 

ds 2 = -N 2 c 2 dt 2 + g i:j (dx l + N l dt) (dx j + N j dt) , 

(i, j = 1,2,3), (1.2) 

which require coordinate transformations be only of the 
types, 

*-►/(*), z'-^fox), (1.3) 
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that is, space-dependent time reparameterizations are no 
longer allowed, although spatial diffeomorphisms are still 
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a symmetry. Then, it is natural, but not necessary, to 
restrict the lapse function N to be space-independent, 
while the shift vector N z and the 3-dimensional metric 
gij in general depend on both time and space, 

N = N(t), N l = N\t,x), ga = g ii {t,x). (1.4) 

This is the project ability condition, and clearly is pre- 
served by the foliation-preserving diffeomorphisms (jl.3l) . 
However, due to these restricted diffeomorphisms, one 
more degree of freedom appears in the gravitational sec- 
tor - a spin-0 graviton. This is potentially dangerous, and 
needs to be highly suppressed in the IR regime, in order 
to be consistent with observations. Similar problems also 
raise in other modified theories, such as massive gravity 

Under the rescaling (jTTTJ) , the dynamical variables 
N, N' 1 and scale as, 



and matter, 



N -> N, N l ->■ e~ 2 N l 



(1.5) 



9ij ' 9ij • 



Note that in [3, [llf , the constant c in the metric (|1.2[) was 
absorbed into N, so that there the lapse function scaled 

as r 2 . 

So far most of the work on the HL theory has aban- 
doned the pro ject ability condition but kept the detailed 
balance [J, |6j, Il2| - |l4j |. One of the main reasons is that 
the detailed balance condition leads to a very simple ac- 
tion, and the resulted theory is much easier to deal with, 
while abandoning projectability condition gives rise to lo- 
cal rather than global Hamiltonian constraint and energy 
conservation. However, with detailed balance a scalar 
field is not UV stable @, and gravitational perturba- 
tions in the scalar section have ghosts 0] and are not 
stable for any given value of the dynamical coupling con- 
stant A [l6j]. In addition, detailed balance also requires 
a non-zero (negative) cosmological constant, breaks the 
parity in the purely gravitational sector [17] . and makes 
the perturbations not scale-invariant 18]. Breaking the 
projectability condition, on the other hand, can cause 
strong couplings [l9| and gives rise to an inconsistency 
theory [H. 

To resolve these problems, various modifications have 
been proposed [21 L In particular, Bias, Pujolas and 
Sibiryakov (BPS) 22| showed that the strong coupling 
problem can be solved without projectability condition 
(in which the lapse function becomes dependent on both 
t and x l ), when terms constructed from the 3- vector 



djN 
N 1 



(1.6) 



are included. Contrary claims can be found in [23j. In 
addition, it is not clear how the inconsistency problem 
[20j is resolved in such a generalization. 

On the other hand, Sotiriou, Visser and Weifurtner 
(SVW) formulated the most general HL theory with pro- 
jectability but without detailed balance conditions [l7j . 
The total action consists of three parts, kinetic, potential 



S = ( 2 J dtd 3 xN^g (L K -C v + C 2 £m) , (1.7) 



where g = det , and 



C K 



KijKV - (1 - K 1 



2 A - R + — (g 2 R 2 + ;/.;/>'.., />'" ) 

+ ^ (/y;/. ,:i + //:,/,' RijRV + fluirjil li: 



+ - [g 7 RV 2 R + g s (V^) (V l i?^)] . (1.8) 

Here C, 2 — 1/WirG, and c denotes the speed of light. 
In the "physical" units, one can set c = 1 \Vj\ . The 
covariant derivatives and Ricci and Riemann terms are 
all constructed from the three- metric gij , while Kij is the 
extrinsic curvature, 

Kij = ^ + ViNj + VjNi) , (1.9) 

where Ni — gijN^ . The constants £, gi (7 = 2, . . . 8) are 
coupling constants, and A is the cosmological constant. 
In the IR limit, all the high order curvature terms (with 
coefficients gi) drop out, and the total action reduces 
when £ = to the Einstein-Hilbert action. 

The SVW generalization seems to have the potential to 
solve the above mentioned problems [24{ , although it was 
found that gravitational scalar perturbations either have 
ghosts (0 < £ < 2/3) or are not stable (£ < 0) 0,111]. 
In order to avoid ghost instability, one needs to assume 
£ < 0. Then, the sound speed c 2 = £/(2 — 3£) becomes 
imaginary, which leads to an IR instability. Izumi and 
Mukohyama showed that this type of instability does not 
show up if \c^\ is less than a critical value [26j . 

It is fair to say, in order to have a viable HL theory, 
much work needs to be done, and various aspects of the 
theory ought to be explored, including the renormaliza- 
tion group flows |27| , Vainshtein mechanism [l5|, [28| , so- 
lar system tests (29|], Lorentz violations [i3], and its ap- 
plications to cosmology 0, [ll| . 

In this paper, we shall study another important issue 
in the HL theory - the problem of singularities, which is 
closely related to the issue of black holes in this the- 
ory [12]. Although we are initially interested in the 
case with projectability condition, our conclusions can 
be equally applied to the HL theory without projectabil- 
ity condition. The extrinsic curvature Kij and the 3- 
dimensional Riemann tensor Rj k i are not tensors under 
the 4-dimensional Lorentz transformations, 



(1.10) 



As a result, in GR one usually does not use them to 
construct gauge-invariant quantities. However, in the 
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HL theory, due to the restricted diffeomorphisms, these 
quantities become tensors, and can be easily used to con- 
struct various scalars. If any of such scalars is singular, 
such a singularity cannot be limited by the restricted co- 
ordinate transformations (|1.3p . Then, we may say that 
the spacetime is singular. It is exactly in this vein that 
we study singularities in the HL theory. In particular, 
we first generalize the definitions of scalar, non-scalar 
and coordinate singularities in GR to the HL theory in 
Sec. II, and then in Sec. Ill we study two representa- 
tive families of spherical static solutions of the HL theory, 
and identify scalar curvature singularities using the three 
quantities K, KyK** and R. In Sec. IV, we present 
our main conclusions and remarks. There is also an Ap- 
pendix, in which we show explicitly that the second class 
of the LMP solutions written in the ADM frame with pro- 
jectability condition in general satisfy the HL equations 
coupled with an anisotropic fluid with heat flow. 

Before proceeding further, we would like to note that 
black holes in GR for asymptotically-flat spacetimes are 
well-defined [3l[ . However, how to generalize such defini- 
tions to more general spacetimes is still an open question 
[32l I33I ] . The problem in the HL theory becomes more 
complicated [IB, 113]) partially because of the fact that 
particles in the HL theory can have non-standard disper- 
sion relations, and therefore no uniform maximal speed 
exists. As a result, the notion of a horizon is observer- 
depedent. 



II. SINGULARITIES IN HL THEORY 

In GR, there are powerful Hawking-Penrose theorems 
31], from which one can see that spacetimes with quite 
"physically reasonable" conditions are singular. Al- 
though the theorems did not tell the nature of the sin- 
gularities, Penrose's cosmic censorship conjecture states 
that those formed from gravitational collapse in a "phys- 
ically reasonable" situation are always covered by hori- 
zons [351 ] . 

To study further the nature of singularities in GR, El- 
lis and Schmidt divided them into two different kinds, 
spacetime curvature singularities and coordinate singu- 
larities |36j ]. The former is real and cannot be made dis- 
appear by any Lorentz transformations (|1.10[) . while the 
latter is coordinate-depedent, and can be made disappear 
by proper Lorentz transformations. Spacetime curva- 
ture singularities are further divided into two sub-classes, 
scalar curvature singularities and non-scalar curvature 
singularities. If any of the scalars constructed from the 
4-dimensional Riemann tensor RZ V \ and its derivatives 
is singular, then the spacetime is said singular, and the 
corresponding singularity is a scalar one. If none of these 
scalars is singular, spacetimes can be still singular. In 
particular, tidal forces and/or distortions (which are the 
double integrals of the tidal forces), experienced by an 
observer, may become infinitely large [371 ] . These kinds 
of singularities are usually referred to as non-scalar cur- 



vature singularities. 

To generalize these definitions to the HL theory, as 
mentioned in the Introduction, both the extrinsic curva- 
ture Kij and the 3-dimensional Riemann tensor can 
be used to construct gauge-invariant quantities, as now 
they are all tensors under the restricted transformations 
(|1.3p . Then, we can see that now there are three kinds of 
scalars: one is constructed from and its derivatives; 
one is from the 3-dimensional Riemann tensor Rjm and 
its derivatives; and the other is the mixture of Kij, Rjki 
and their derivatives. Therefore, we may define a scalar 
curvature singularity in the HL theory as the one where 
any of the scalars constructed from K%j, Rl kl and their 
derivatives is singular. A non-scalar curvature singular- 
ity is the one where none of these scalars is singular, but 
some other physical quantities, such as tidal forces and 
distortions experienced by observers, become unbounded. 
Coordinate singularities are the ones that can be limited 
by the restricted coordinate transformations (|1.3|) . 

Several comments now are in order. The scalars con- 
structed from the 3-dimensional tensors Kij and R^i and 
their derivatives include all scalars constructed from the 
4-dimensional RZ V \ and its derivatives. Thus, according 
to the above definitions, all scalar singularities under the 
general Lorentz transformations (II . 10j> are also scalar sin- 
gularities under the restricted transformations (|1.3j) . but 
not the other way around. In this sense, scalar singulari- 
ties in the HL theory are more general than those in GR. 
One simple example is the anti-de Sitter Schwarzschild 
solutions, which are also solutions of the SVW general- 
ization with £ = 0, as in this case the 3-dimensional Ricci 
tensor Rij vanishes identically, and the contributions of 
high order derivatives of curvature to the potential Cy 
are zero, as can be seen from Eq. (|1.8[) . However, as 
shown in the next section, the corresponding two scalars 
K and KijK ij all become singular at r = (3A//| AI) 1 / 3 . 
This singularity is absent in GR [3l| . 

In 3-dimensional space, the Weyl tensor vanishes iden- 
tically, and the Riemann tensor is determined alge- 
braically by the curvature scalar and the Ricci tensor: 

Rijki — gikRji + gjiRik — gjkRu — guRjk 

-\{9ik9ji~ gu9jk)R- (2.1) 

Therefore, the singular behavior of the scalars made of 
the 3-dimensional Riemann tensor i?* fei may well be rep- 
resented by the 3-dimensional curvature scalar R. 

III. SINGULARITIES IN SPHERICAL STATIC 
SPACETIMES 

The metric of general spherically symmetric static 
spacetimes that preserve the ADM form of Eq. (|1.2j) 
with the projectability condition can be cast in the form 

E3, 

ds 2 = -dt 2 + e 2l/ (dr + e^-^dt) 2 + r 2 dn 2 , (3.1) 
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where /i = /i(r), v = v(r). Then, we have 



R = 
K = 



r- 



2rV - (1 - e 2v ) 



(3.2) 



where v 1 = dv j 'dr, etc. It is interesting to note that for 
the metric (13.11) we have 



(3.3) 



dj = = e iik RuVjB! k , 
where Ca is the Cotton tensor, defined as 



4 1 



(3.4) 



As a result, the HL theory with detailed balance pj 



S. 



HLd 



dtdx 3 N v Kjl — 



K l3 K ij - \K 2 



«V (A W R - 3A 2 W ) K 2 /i 2 (1 - 4A)i? 2 



(1-3A) 



32(1 - 3A) 



A. ,2 



2^° ijG 



(3.5) 



can be effectively considered as a particular case of the 
general SVW action (|1.7I) with 



G = 



327TC 2 

4A - 1 



kV 2 A^ 
16(1 -3A)' 



A = 



3A 



tr 



'<% 93 



1 - 3A 



4A W A w 
£ = 1 - A, 34 = 05 = •■• = 5s = 0. 



(3.6) 



It should be noted that these relations are valid only 
when the conditions of Eq. (|3.3p hold. In general, these 
two terms do not vanish and violate parity, while the 
SVW action always preserves it. It must not be confused 
with the parameter fi used in the action (|3.5p and the 
metric coefficient used in (13.11) . 

To study singularities in the HL theory, in the rest of 
this paper we shall restrict ourselves to two representative 
cases, the (anti-) de Sitter Schwarzschild solutions and 
the solutions found by Lu, Mei and Pope (LMP) [6J. 



A. (Anti-) de Sitter Schwarzschild Solutions 

The (anti-) de Sitter Schwarzschild solutions are given 

by m 

»=Hv + r 2 )> » =0 - (3 - 7) 



When A > 0, it represents the de Sitter Schwarzschild 
solutions, and when A < it represents the anti-de Sitter 
Schwarzschild solutions. As mentioned previously, they 
are also solutions of the SVW generalization with £ = 0. 
Inserting the above into Eq. (|3.2l) . we find that 



R 
K 



= 0, 



^3Af 



Ar 



12r 3 
3M + Ar 3 



1/2 



3M - 2Ar 3 
3M + Ar 3 



12r 3 



/3M - 2Ar 3 
V 3M + Ar 3 



(3.8) 



Clearly, when A > 0, K and KijK 1 ^ are singular only at 
the center r = 0. However, when A < 0, they are also 
singular at r = ta = (3M/IAQ 1 / 3 . In contrast to GR, 
this singularity is a scalar one, and cannot be removed 
by any coordinate transformations given by Eq. (|1.3j) . 

In GR, the (anti-) de Sitter Schwarzschild solutions are 
usually given in the orthogonal gauge, 



with 



ds 2 



4- = 



. e a*(r) dT a + e a»(r) <ir a + , 



■ 2 dtt 2 , 



-<E> = - In 
2 



2M 
r 



-Ar 



(3.9) 



(3.10) 



Clearly, the metric (|3.9[) does not satisfy the pro- 
jectability condition, and its coefficient g rr is singular 
at e 2 ®( r EH) — o But, in GR this is a coordinate singu- 
larity, and all scalars made of the 4-dimensional Riemann 
tensor and its derivatives are finite at r — teh- 

It is interesting to note that in the orthogonal gauge 
p.9[) . K and K^K % ^ all vanish, as can be seen from Eqs. 
(|1.9[) . while the 3-dimensional Ricci scalar is given by 
Rorth — — 2A, where R or th denotes the quantity calcu- 
lated in the orthogonal gauge (|3.9[) . 

In (38|, it was showed explicitly that metric (I3.9[) is 
related to metric (13. 1[) by the coordinate transformations, 



t- 



J" Ve- 2 * - 1 



"dr, 



(3.11) 



under which we have 



$(r) = v{r) - -ln^l- e 2 ^, 



*(r) = |ln(l-e 2M 



(3.12) 



or inversely, 
1 



H= -ln(l-e^), zy = $(r) + *(r). (3.13) 



Note that the coordinate transformations p. lip are 
not allowed by the foliation-preserving diffcomorphisms 
(ll.3p . In addition, since K, K^K^ and R are not scalars 
under these transformations, it explains why we have 
completely different physical interpretations in the two 
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different gauges, defined, respectiveiy, by Eqs. (|3.f| and 
(|3.11[) . In fact, in the framework of the HL theory the two 
different gauges represent two different theories - metric 
(|3.1[) represents a HL theory with projectability condi- 
tion, while metric (|3.11l) represents a HL theory without 
projectability condition. 



B. The LMP Solutions 

The LMP solutions were originally found for the HL 
theory with detailed balance but without projectability 
conditions in the form (|3.9j) . There are two classes of 
solutions, given, respectively, by 



\n(l + x 2 ), 



(3.14) 



for any ^(r), and 
1 



$ = -|ln(l + :r 2 - ax a± 



* = ~P± ln(ar) + | In ( 1 + x 2 - ax a± ) , (3.15) 



where x = \J— Aw r, and 



2A - V6A - 2 
A - 1 



2(2A- 1) 



2A + V6A-2 
A - 1 ' 



2A + V6A - 2 

/3± = 2a± - 1, 



(3.16) 



with a being an arbitrary real constant. Fig. [T]schemati- 
cally shows the curves of a± vs A, which shows that there 
is no solution for the a+ branch when A = 1. 

Solutions given by Eqs. (|3.14j) and (|3.15j) shall be re- 
ferred to, respectively, as Class A and B solutions. Par- 
ticular values of a± are 



a_(A) 



a+(A) 



f-1, 


A = 


1/3, 


o, 


A = 


1/2, 


< V2, 


A = 


1, 


2/3, 


A~ 


1.375 


1, 


A = 


3, 


I 2, 


A = 


00, 




A = 


= 1/3, 


1 -00, 


A = 


= 1-°, 




A = 


= l +0 , 


I 2, 


A = 


= 00, 



(3.17) 



which are useful in the following discussions. 

In the orthogonal gauge (|3.9p , as mentioned previously, 
the extrinsic curvature Kij vanishes identically, while the 
3-dimensional curvature for the above two classes of so- 
lutions can be written as 



R 



■th 



-{a( 



1 + a±)x c 



3x 2 }. 



(3.18) 



When a = 0, the corresponding R or th is for the solution 
(|3.14[) , which reduces to a constant. That is, in this case 




FIG. 1: The functions a±(A) defined by Eq. ((3TT6J. 



all the three scalars K, K^K'^ and R are finite. When 
a/0, it is for the solution (|3.15[> . which is singular only 
when a± < 2 at r = 0. From FigfT] we can see that 
a_ is always less than two for a finite A. Therefore, this 
branch of solutions is always singular at the center. a + 
is always less than two for 1/3 < A < 1 and greater 
than two for A > 1. Therefore, R or th is finite for the ab- 
solutions with A > 1 for any r, including the center r = 0, 
while it is singular at the center for the a+ solutions with 
1/3 < A < 1. 



1, Class A Solutions 

Transforming the above solutions into the canonical 
ADM form (j3~Tj) . we find that for the solution (|3~T4]) , we 
have 



fi = — oo, v = — — In ^1 — A\yr 2 



(3.19) 



where in writing the above expressions, we had chosen 
"J = 0. Then, the corresponding three scalars are given 
by 



K = K i3 K l] = 0, 



i? = 6A 



w ■ 



(3.20) 



which are all finite. As shown in [38|, this solution is also 
a vacuum solution of the SVW generalization (ll.7[) with 
a non- vanishing constant curvature k — 6Avi/. Since it 
does not depend explicitly on the coupling constant £, it 
is a vacuum solution for any given £, including £ = 0. 
For detail, we refer readers to Sec. V of (38j . 
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2. Class B Solutions 



For Class B solutions (|3. 15p , they can be written in the 
canonical ADM form (13.11) with 



/' = 

where 

A = 1 



-In A, v 
2 



p 2(l-2a±) 



;i - 2a±) \n{x), 



, a .4(l-a ±)+aa .2-3a ±> 



(3-21) 



(3.22) 



Clearly, to have real solutions, we must assume A > 0. 
It should be noted that unlike Class A solutions, this 
class of solutions do not satisfy the vacuum equations of 
the HL theory with projectability condition, due to the 
fact that the HL actions are not invariant under the co- 
ordinate transformations (|3 . 1 1 [) . As shown in Appendix, 
they can be interpreted as solutions of the HL theory 
with projectability condition coupled with a spherical 
anisotropic fluid with heat flow. The properties of sin- 
gularities of these quantities can be well represented by 
the three scalars K, K^K 1 ^ and R, as can be seen from 
Appendix, so in the following we shall not consider them 
specifically. 

Inserting the above solutions into Eq. (I3.2p we find 
that 



R 
K 



A{i-(i + 2/?± )^ ± } 



2A 1 /2 a; 3-2 

( 8 A - X 2 S 2 ) 



K W A- x^S) 
2x2 ^ (3.23) 



where 



5 = 2(l-2a±)a; 1 - 4Q ± +4(l-a±)a; 3 - 4Q± 

- a(2-3a ± )x 1 - 3a± . (3.24) 

To study these solutions further, it is found convenient 
to consider the two branches a± separately. We shall use 
P = P± to denote the a± branches. 

Case i) f3 = f3 + , 1/3 < A < 1: In this case we have 
a + < — 1 and j3 + < —3. Then, from the expression (I3.22[) 
we find that 



x = 0, 
x — » oo. 



(3.25) 



Thus, there must exist a point x = x s at which we have 
A(x s ) — 0. Since A > 0, we must restrict the solutions 
to the range < x < x s (or equivalcntly, < r < r s , 
where r s = x s j \J— ~Kw-)- Fig. [2] shows this case. From 
the expressions (|3.23p we find that all these three scalars 



diverge at the center 



0, while K and KuK %3 di- 



verge at r = r s . Therefore, in the present case there are 
two scalar curvature singularities, located, respectively, 
at r = and r = r s . 



^ A « cc>0 




FIG. 2: The function A(x) defined by Eq. (ET22")| for the a+ 
branch with 1/3 < A < 1. All three scalars R, K and KijK 13 , 
defined by Eq. (|3.23[) . become unbounded at r — 0, while only 
K and KijK' J diverge at r = r s . 



Case ii) j3 = f3+, A > 1: In this case we have a+ > 2 
and /3+ > 3, where the equality holds only when A = oo. 
Then, from the expression (|3.22p we find that 



A = 



r-oc 



x 
x 



0, 

> oo. 



(3.26) 



FigEl shows this case, from which we find that for any 
given A > 1, there always exists a minimum r s > so 
that A(r) > for r > r s , where r s is the solution of 
A(r) = 0. Then, from the expressions (|3 . 23[) we find 



that K and K^K 13 diverge at r 



while R diverges 



as r — >• oo. Therefore, in the present case there are also 
two scalar curvature singularities, located, respectively, 
at r = r s and r — oo. 

Case iii) (3 = 1/3 < A < 1: In this case we have 
— 1 < a_ < 1/2 and —3 < {3- < 0, where the equality 
holds only for A = 1/3. Then, from the expression (I3.22[) 
we find that 



A = 



I -ex 



X 
X 



0, 

► oo. 



(3.27) 



from Figf4l we find that for any given A in this range, 
there always exists a maximum r s > so that A(r) > 
for < r < r s . Then, From the expressions (|3.23[) we 
find that all three scalars R, K and K^K 13 become 
unbounded at the center, while only K and K^K 13 di- 
verge at r = r s . Therefore, in the present case there are 
two scalar curvature singularities, located, respectively, 
at r — and r = r s . 

Case iv) /3 = /3_ = 0, A = 1: In this case we have 
a_ = 1/2, and 



ax 



1/2 



x 2 . 



(3.28) 



Thus, to have A non-negative, we must assume a > 0. 
Then, A > for < r < r s = a 2 ' 3 , where A(x) = at 
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H Mr) a > a c 




r 



FIG. 3: The function A(x) denned by Eq. (|X22"j) for the a+ 
branch with A > 1, where a c is the value of a, for which 
A (a:) > for a < a c . K and KijK 13 diverge at r = r s , while 
R diverges as r — > oo. 



i i A W a > a c 




FIG. 4: The function A(x) defined by Eq. (f3722]l for the a_ 
branch with 1/3 < A < 1, where a c is the value of a, for 
which A'(x) < for a < a c . i?, K and KijK 13 diverge at the 
center r = 0, while only and KijK lJ diverges at r = r s . 

both r = and r = r s . At the center, all three scalars 
R, K and KijK^ become unbounded, while only K and 
KijK t; > diverge at r — r s . 

Case v) /3 = /?_, 1 < A < 3: In this case we have 
1/2 < a_ < 1, < /3_ < 1, and 

A = (~°°' 31 = °' (3.29) 
I — oo, x — > oo. v ' 

Fig JS] shows the general properties of A(r), from which 
we can see that for any given A, there always exists a crit- 
ical value a c , for which A(r) > is possible only when 
a > a c . In the latter case, A(r) = always has two 
positive roots, say, T\ and r-i- Without loss of general- 
ity, we assume r<i > r% > 0. When n < r < ?'2, A(r) 



il A(r) 



a > a c 




FIG. 5: The function A(z) defined by Eq. (1^22)) for the 
a_ branch with 1 < A < 3, where a c is the value of a, for 
which A(x) > is possible only when a > a c . K and KijK 13 
diverge at the two points r = ri,2, while i? remains finite 
there. 

is non-negative. At the two points r = n and T2 we 
have A(r,*) = 0, and Eq. p.23p shows that both K and 
KijK 1 ^ become unbounded at these points, while R re- 
mains finite. Therefore, in the present case there are also 
two scalar curvature singularities, located, respectively, 
at r = ri and r = T2 for a > a c . Solutions with a < a c 
are not physical. 

Case vi) (3 = j3- = 1, A = 3: In this case we have 
a_ = p— = 1, and 

A=-- — = ( _0 °' x = ' (3 30) 

Clearly, to have A > 0, we must assume a > 0. Then, 
for r>r s = l/(^/—Aw a), A is non- negative [See Fig|6l 
Curve (a)]. At r = r s , A(r) = 0, and Eq. ([3T23]) 
shows that both K and KijK^ become unbounded at 
this point, while R remains finite. As r — > oo, A(r) — > 0, 
and K and K^K^ become unbounded again, while R 
still remains finite. Therefore, in the present case there 
are also two scalar curvature singularities, located, re- 
spectively, at r = r s and r = oo. 

Case vii) {3 = /3_, A > 3: In this case we have a- > 
1, (3- > 1, and 

A = {-~' X = °' (3.31) 
I +1, x —> oo. v ; 

FigEfb) shows the general properties of A(r), from which 
we can see that for any given A, there always exists a 
point r — r s where A(r s ) = 0. When r > r s , A is 
always positive. Eq. (|3.23p shows that now both K and 
KijK^ diverge at r — r s , while all three scalars become 
unbounded as r — > oo. Thus, in the present case there are 
two scalar curvature singularities, located, respectively, 
at r = r s and r = oo. 
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A A(r) (b) 




FIG. 6: The function A(x) denned by Eq. (f3~22j) for the q_ 
branch with (a) A = 3, a > 0; and (b) A > 3 for any a. 



IV. CONCLUSIONS 

In this paper, we have studied singularities in the HL 
theory, and classified them into three different kinds, the 
scalar, non-scalar, and coordinate singularities, following 
the classification given in GR [36j . Due to the restricted 
diffcomorphisms (|1.3I) . the number of the scalars that 
can be constructed from the extrinsic curvature tensor 
Kij, the 3-dimensional Riemann tensor and their 
derivatives is much larger than that constructed from 
the 4-dimensional Riemann tensor RZ V \ and its deriva- 
tives. The latter is invariant under the general Lorentz 
transformations (|1.10[) . As a result, even for the same 
spacetime, it may be singular in the HL theory, but not 
singular in GR. One simple example is the anti-de Sitter 
Schwarzschild solution written in the ADM form (|3.ip . 
This solution is a solution of both HL theory and GR. 
However, in the HL theory, there are two scalar singu- 
larities located, respectively, at the origin r = and 
r = (3A//IAI) 1 / 3 , as the two scalars K and /\,,/\ ■' be- 
come singular at these points. It is well-known that the 
scalar singularity at r = (3M/IAI) 1 / 3 is absent in GR. 
This is because both K and KijK^ are not scalars un- 
der the general Lorentz transformations (jl.101) . Thus, 
even they are singular at this point, it only represents a 
coordinate singularity, and all scalars constructed from 
the 4-dimensional Riemann tensor and its derivatives are 
finite. On the other hand, due to the restricted transfor- 
mations (jl.3[) , K and KijK'^ are scalars in the HL the- 
ory, and once they are singular, the resulting singularity 
cannot be transferred away by the restricted transforma- 
tions. As a result, in the HL theory it represents a real 
spacetime singularity. 

With the above in mind, we have studied the LMP so- 
lutions Q, and found that their singularity behavior in 
the orthogonal frame defined by (|3.9|) is different from 
that in the ADM frame defined by (|3.ip for the second 



class of the LMP solutions. In particular, in the orthog- 
onal frame, vanishes, so do K and KijK^ , while the 
3-dimensional Ricci scalar R [cf. Eq. (|3.18[) ] can be sin- 
gular at the origin or infinity, depending on the choice of 
the parameter A. However, in the ADM frame at least 
one of the three scalars K, K l i and R is always sin- 
gular at two different points, either r = and r — r s > 0, 
or r = r\ and r = T2 with ri > r\ > 0, or r = r s > and 
r = oo, depending on the choice of the free parameter A, 
where r s is a finite non-zero positive constant. This dif- 
ferent singular behavior originates from the fact that the 
two frames are related by the coordinate transformations 
p. lip , which is not allowed by the foliation-preserving 
diffeomorphisms (11.31) . or in other words, K, KijK^ and 
R are not scalars under such transformations. In fact, 
in the framework of the HL theory, the two frames actu- 
ally represent two different HL theories, one is with the 
projectability condition, while the other is without. In 
particular, the second class of the LMP solutions in the 
orthogonal frame satisfy the vacuum HL equations, while 
in the ADM frame they satisfy the HL equations coupled 
with an anisotropic fluid with heat flow, as shown explic- 
itly in the Appendix. 

Our above results show clearly that the problem of 
singularities in the HL theory is a very delicate prob- 
lem, due to the restricted diffeomorphisms (ll.3p . which 
preserve the ADM foliations (jl.2[) . Further investiga- 
tions are needed, in particular, in terms of the strength 
of these singularities. In the examples studied in this 
paper, all singularities indicated by the two scalars K 
and KijK 1 ^ at r = r s > 0, including that of the anti-de 
Sitter Schwarzschild solution, seems weak in the sense 
of tidal forces and distortions experiencing by observers. 
Therefore, it is not clear whether or not the spacetime is 
extendable across such a singularity 37]. 

Finally, we would like to note that the ADM form (|3.1[) 
can be considered as a particular case of the HL theory 
without projectability condition. Therefore, restricting 
ourselves only to the HL theory without projectability 
condition does not solve the singularity problem occur- 
ring at r = r s . 
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Appendix: The Generalized LMP Solutions with 
Projectability Condition 

The second class of the LMP solutions with pro- 
jectability condition takes the form of Eq. (13.11) with 
H and v given by Eq. p. 211) . Unlike the first class, 
this one does not satisfy the HL vacuum equations with 
projectability condition, due to the fact that HL actions 
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are not invariant under the coordinate transformations 
(|3.11[) . In particular, under these transformations, we 
have 



(A.l) 



For spherically symmetric static solutions, the extra term 
SRij in g eneral gives rise to an anisotropic fluid with heat 
flow [38| ], possibly subjected to some energy conditions 
[3lj |. In particular, for the LMP solutions of Eqs. fl3.1[) 
and (|3.21l) . the energy density J*, defined by 



(A.2) 



is given by the Hamiltonian equation, 

J d 3 x^/g(£ K + C V ) =8nG J d 3 x y ^J t : (A.3) 



where 



\^x 2 A' 2 - 8(1 - £)xAA' 



C v = 



x 2 



Ax 2 A 
l-20xA 2 }, 

2 + 3a; 2 + 2(l - Aa±)x 2f}± 



A 



ir 



■{[8^o4-8(l-0(l-a±)+3] 



x 8(l-a±) 

+2 [4(1 - 0(1 ~ o±) - l]a; 2(1 ~ 2a±) 
+ (1-20x 4 ( 1 - 2q± '}, 



(A.4) 



where A is given by Eq. ([3T22"]) . and A' = dA/dx. The 
quantity u, defined by 

J* = -JV^=e-^+")(«,0,0), (A.5) 



which is related to heat flow [38j, is given by 
A w e 2 (^> ' 



32ttGx 2 A 2 



ix 



2a;AA" + a;A' - 2/3±AA' 



4A £xA' - 2(1 - 0/?±A - 8£A 



(A.6) 



On the other hand, the corresponding stress part T iJ 
defined by, 



2 S(y/gC M ) 



(A.7) 



can be written in the form, 



where fly = 5$6$ + sin 2 95? Sf, and 



(A.. 



647rGa; 2 A 2 



4xAA" - 3a; A' + 4/3± A A' 



+8A 



xA' - 2(1 - 0/3±A + 8(1-4£)A 



~8^G 
167rGa; 2 A 2 



(1 - 0x 2 AA" + -£a; 2 A 



2 A' 2 



-(1 - 0/?±a;AA' + 2(1 - 2£)(xA' + f3±A) A 



A, 



87rGa; 2 ' 



(A.9) 



with 



Aw 

2a; 2 



2 - 2x- 2 ^ - toPV-M 
Aw 



(19 - 32^) 



4a ,4-2/3± 

-4(3- 50(2 -30 ± )0± 

-2(7 -U^x' 2 ^ - (5-8^0;^ 



2 P± 2a; 2 - 4 ' 3 ± 



(1 + 20 



+2(ll-160^±-2^7-6^ ± )^ 



-2(l-/3±)x- 2/3± +(l-2£)a 



-4/S± 



.(A.10) 



Clearly, to interpret the above as an anisotropic fluid 
with heat flow, some energy conditions [3l[ might need to 
be imposed. Since in this paper we are mainly concerned 
with the nature of singularities, we shall not study these 
conditions here. 
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